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ABSTRACT 


This paper treats, from a statistical mechanics view- 
point, the justification for acknowledging satellite posi- 
tional uncertainty through in -track error buildup. In ad- 
dition, by relating the angular momentum to a quantity 
defined as the virial parameter, through a derived quan- 
tity defined as the Vector Density function, it is shown 
that for any conservative system, the periodic variations, 
and hence the cross track error, will time average to 
zero. This is a direct consequence of the conservation 
of angular momentum. 



SUMMARY 


The present paper treats, from a statistical mechanics viewpoint, 
the justification for acknowledgment of satellite positional uncertainty 
through in -track error buildup. The equations of motion of the Hamilton - 
Jacobi theory specify a region or volume in phase space in which the 
system point moves. Since the system is in statistical equilibrium, 
it will initially have a specified energy, and the theorem of the integral 
invariants of Poincare guarantees that volume in phase space in which 
the system point moves with an initially specified energy will be con- 
stant in time, although forces influencing the system will cause a dis- 
tortion or elongation of the volume with time. Therefore, in traversing 
this distorted volume, an error factor in the equations of transforma- 
tion will necessitate a change in the constant of the motion which relates 
to the time of perigee passage and, hence, the mean anomaly. Because 
the secular part of the mean anomaly is by far the largest portion, any 
error accumulated here can be brought out in terms of an in-track pro- 
jection and will be mainly secular in nature. In addition, by relating 
the angular momentum to a quantity defined as the virial parameter, 
through a derived quantity defined as the Vector Density Function, it is 
shown that for any conservative system the periodic variations (of the 
orbit), and hence the cross-track error, will time average to zero. 
This is a direct consequence of the conservation of angular momentum . 
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by 
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INTRODUCTION 

If we consider a 2n-dimensional Cartesian space formed of coordinates q x • • • q n , p : • • • p n , 
and known as phase space, the complete dynamical specification of a mechanical system will be 
given by a point in such a space. The values of q and p at any time t can be obtained from their 
initial values by a canonical transformation that is a continuous function of time; i.e., there exists 
a canonical transformation from the values of the coordinates and momenta at any time t to their 
initial values. Obtaining this transformation is equivalent to solving the problem of the system 
motion 0 While the exact motion of any system is completely determined in classical mechanics by 
the initial conditions, in practice these conditions are often only incompletely known. In addition, 
the presence of nonconservative forces and uncertainties in fundamental data introduces further 
errors into the solution for the motion of such a system. By introducing certain techniques of 
statistical mechanics, one can make predictions about certain average properties by examining 
the motion of a large number of identical systems. The values of the desired quantities are com- 
puted by forming averages over all systems in the ensemble. Statistical techniques can be justi- 
fied by showing that the time average of any quantity pertaining to any single system of interest 
would actually agree with the ensemble average for that quantity calculated for all members of 
the corresponding ensemble. The postulate leading to this conclusion is the ergodic hypothesis 
that states that the phase point for any isolated system will pass in succession through every 
point before finally returning to its original position in phase space. 


STATEMENT OF THE PROBLEM 

Consider the motion of the system point in phase space for the case of an artificial earth 
satellite. For such points, Liouville's Theorem states that the density of systems in the neighbor 
hood of some given system in phase space remains constant in time. That is, 

dD _ <9D ^ ( d D . _dD . \ 

dt dt + / , ^dp 4 p i + dqi q i/ ’ 

i = 1 


1 


or, 


dD 

dt* 


dV 

' dt + 


[D,H] 



where D is the density, H is the Hamiltonian, Nis the number of particles, and Pi and q. are the 
generalized or canonical momenta and position coordinates, respectively. These are given by 


and 


r i (9H _ . 

= "TqT = Pi • 


r 1 dti 

= = q i ’ 


When the ensemble of systems is in statistical equilibrium, the number of systems in a given 
state must be constant in time. Then [d,h] = 0, or the Poisson bracket vanishes. Therefore, we 
have equilibrium by choosing the density to be a function of the constants of the motion of the 
system. For conservative systems, D can be any function of the energy. The density is a constant 
equal to zero for the system with a given energy considered here. 

Now canonical transformations are defined as having the property of preserving the form of 
Hamilton’s equations of motion under transformation. There are other expressions invariant 
under canonical transformation* The theorem of Poincare states that the integral 

Jl = jj2Z dq ' dp - 


is invariant under canonical transformation. The proof of invariance can be reduced to showing 
that the sum of the Jacobians is invariant (Reference 1). In a similar fashion, the integral invari- 
ants can be extended to 


J 


n 



■■■ dq n d P j • • • dp n , 


where the integral can be extended over an arbitrary region in phase space. The invariance of J n 
is equivalent to the statement that volume in phase space is invariant under canonical transforma- 
tion. From this it follows that the volume in phase space is constant in time. Now the motion of a 
system point is simply a particular contact transformation of the canonical coordinates. Conse- 
quently the change in time of a region of phase space can be represented by a contact or canonical 
transformation; that is, the contact transformation determines motion within the region of constant 
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volume. From the Hamilton -Jacobi theory, the motion of an artificial earth satellite for a con- 
servative system can be described by 


‘ = 

^ = 

and 

^ -- 

where a l9 a 2 , and a 3 are related to the total energy, total angular momentum, and polar component 
of angular momentum, respectively; fi 1 , ft 2 , and are related to the time of perigee passage, 
argument of perigee, and geocentric right ascension, respectively; and w is Hamilton's character- 
istic function (References 2 and 3). 

These equations of motion, then, specify a region or volume in phase space in which the sys- 
tem point moves* In the course of time, the volume may take on different shapes, but as long as 
the system remains a conservative one the volume remains constant. If now any of the earth's 
geodetic constants are incorrectly specified or ignored in the representation of the gravitational 
potential 


<9 W 
da x 


J _ W 

d CL 2 


d W 


V 



">(c, 

n = 1 m = 0 


cos mA. + S 


in m\) 


where /z is the product of the gravitational constant and the planet's mass, r e is its equatorial 
radius, r and & are the p lan eto centric distance and latitude of a field point, C nm and S n m are the 
coefficients of potential, \ is the angle of east longitude, and P n m (cos <9) is the associated Legendre 
functions, then the path of the system point in phase space will deviate more and more in time 
from one governed by the exact representation of the potential given above. In addition the initial 
volumes described will differ, and each of these will again evolve separately in time. The system 
point for this orbit which is bounded and not closed will therefore not arrive at the place and time 
in phase space determined by the correct solution* This error is manifest principally in the 
equation of motion 


t + /3 X 


aw 
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and indicates that the time of arrival of the satellite at perigee is changed. With each revolution, 
the difference in the times of perigee passage increases, and the ’"position” of the system point in 
phase space becomes more removed from that of the exact solution. Therefore, there is a stretch- 
ing or elongation of system point positions in phase space, corresponding to a time and position 
error for the satellite. Since the time of perigee passage is related to the mean anomaly by the 
relation M = n (t +/3J, where n is the mean motion, then the error is accumulated directly by the 
mean anomaly, a fundamental orbital element. It can be shown that the secular part of the mean 
anomaly is at a minimum, on the order of J 2 times the periodic part, where J 2 is the coefficient of 
the second zonal harmonic of the earth’s gravitational potential (Reference 2). 


ANGULAR MOMENTUM AND THE VECTOR DENSITY FUNCTION 

Referred to a fixed origin, the angular momentum of a particle of mass m is given by 
L = (7 x p), where the momentum p = mv, and r is the position vector of the mass. For unit mass, 
this is given in component form by 

L = (rxv) = i ( v z - r z v y ) + T(r z v x - r x v z ) + k(r x v y - r y v x ) , 

with respect to a set of Cartesian axes. Let us write the virial parameter G, y = r e v y , where ft 
and y (which equal 1, 2, or 3 exclusively), refer to a particular component x, y, or z. Now let us 
define a Levi-Civita type Vector Density ~% aBy , where a, /3, and 7 are 1, 2, or 3 exclusively, s is 
the unit vector T, 7, or E, according to whether a = 1, 2, or 3, respectively; positive or negative 
for a/jy an even or odd permutation of 123; and zero if any two subscripts are equal. For example, 
for a/3y = 123, 1> 1J3 = +T; if a/3y= 213, ^ 213 = ~T. Using this, the angular momentum becomes 

L -- ?[G y ,-G zy ] + T[G zx -G X2 ] +k[G xy -G yx ] , 

or simply, 

L - G /3y ' 

a,p,y-l 


Therefore, the angular momentum for the system can be expressed as individual sums of the virial 
function through the use of the vector density, § a/3y . 


POSITIONAL UNCERTAINTY 

From the above, it is important to consider how any error in orbital elements, especially the 
very sensitive mean anomaly, is reflected in the satellite's positional uncertainly. Then the 
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in -track position error can be defined as 


AT 


v • A r 


where r = lx + J y + kz and v = ix + jy + kz. Thus the in -track error is the error in position 
given by Ar, proj ected along the unit tangent velocity vector. In addition the error normal to the 
orbit plane is given as 


an 


L • Ar 

|L| 


where 


L - ( r x v) . 

Here the normal or cross-track error is the error in position Ar, projected in the direction of the 
unit angular momentum vector. Also I U is the magnitude of the constant angular momentum of 
the orbit. If now the time average of these quantities is examined, then for the cross -track error 

_1_ f T -A * , _L f T JL L - Ar 

r ] o dt ^ dt - r ] o dt (El dt • 

where r is some time interval. Since Lis a constant of the motion, we can write 



. \ 

3 

\ 

f T d 1 

1 « <L ' ar > d ' ‘ r|E| 

o 

y \ ^ a /3y G /3r 

a. , ft , y- 1 

1 


= 7^ + ([ G :*-G xz ]Ay) + ([Gx,-Gy>)}o T • 

= 7j7y{([ G yz (nA*(r) _ G yz ((»Ax(0)] - [G 2y (r)Ax(r) -G 2y (O)Ax(O)]) 

+ ([G zx (r-)Ay(r) -G 2X (O)Ay(O)] - [G X2 (r)Ay(r) -G X2 (O)Ay(O)]) 

+ ([G xy (t)Az(t) -G xy (O)Az(O)] - [G yx (r)Az(r) -G yx (0)Az(0)])| . 
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Now it can be shown mathematically that Ax, Ay, and Az are functions of the changes in the six 
elements of the orbit, and for an elliptic orbit the variation in these elements will be complex 
periodic functions of the argument of latitude, true anomaly, eccentric anomaly, and inclination, 
and will hence be bounded (see appendix). Therefore, Ax, Ay, and Az will exhibit this behavior and 
also be bounded. This is fairly obvious since the difference in Cartesian components between two 
ellipses will not be infinite. In addition, since the coordinates and velocities of the satellite remain 
finite, there is an upper bound to . Then by choosing r sufficiently long, each of the products, 

T ( G /3y ^ r a) ’ 


where Ar a = Ax , A y , or Az for a = 1, 2, or 3, respectively, and a \ /3 \ y, can be made as small as 
desired. Therefore, in this limit the time average for the error normal to the orbit with an 
angular momentum that is a constant of the motion is said to vanish. 

For the in -track position error, since AT does not involve L or products of the form v but 
rather components of the velocity and the magnitude of velocity, the above arguments do not apply; 
hence, 


_1 

T 


r 1 

Jo 


d 

~dt~ AT dt 


f 0 . 


Therefore, errors of the normal will be periodic in nature with both short and long periodic effects 
present, at times causing distortion of the error curve. The in-track error curve will then be 
secular -like and can have periodic effects added. It should become larger than the cross-track 
error and grow without bound (Reference 4). This result seems reasonable in view of the fact that 
the major portion of the error buildup appears in the time of perigee passage, or mean anomaly, 
whose secular part is much larger than the periodic part. 


REMARKS 

The effects of the various terms of the associated Legendre functions (Reference 4) are briefly 
summarized here. 

1. For the zonals, the effect is either secular or periodic, depending on whether n is even or 
odd. 

2. The sectorial harmonics (n = m) always produce secular -like changes. 

3. For the Tesseral harmonics (n / m), the C n m terms induce a large periodic change with a 
secular -like effect added, and the S n m terms will induce a secular effect with a periodic 
variation added. 
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Now, if nonconservative forces of a quasistatic nature are present, that is, if they draw 
energy very slowly, then the above discussions will apply to a great extent. For example, if a 
small amount of energy is withdrawn from the system only at the point of perigee passage, the 
above results are applicable after the orbital constants have been redetermined. In this way, the 
angular momentum will be a constant of the motion for a given revolution. If the total energy of 
the orbit is being constantly changed, then all time-average arguments do not apply, but the statis- 
tical mechanical considerations are still valid. Consequently, for a nonconservative system, as 
soon as the nonconservative forces are removed, the above time-average argument applies again. 
Since an exact classical description for the effects of these nonconservative forces is not availa- 
ble, then from a statistical mechanical standpoint the discrepancy between the trajectories of the 
phase points of the actual and computed system will be even greater. Presently, therefore, a rea- 
sonable determination of the positional uncertainty of a satellite can be made providing the change 
in the time of perigee passage can be accurately determined. 

In addition, for conservative systems, as a consequence of the conservation of angular momen- 
tum, periodic variations in the orbit, and therefore the cross-track position error, will time 
average to zero. The secular changes and hence the in-track position error will grow steadily. 
This error is reflected principally in the time of perigee passage or the mean anomaly. 

While the purpose of this paper has been to study the physical meaning of effects caused by 
the errors from gravity models, nonconservative forces, etc., on the satellite trajectories, the 
problem now reduces to that of determining, if possible, an accurate estimate or error bound for 
each of these sources and the resulting magnitude of this error bound projected through the mean 
anomaly and in -track position. This will constitute the subject matter for a future investigation. 
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Appendix 

Discussion of Oscillatory Nature of Error in Satellite Coordinates 


It was pointed out in the section on positional uncertainty that the error bound in the position 
of a satellite could be shown to be a complicated oscillatory function of several orbital elements. 
To do this, it should be assumed that any errors or uncertainty in the potential model of the earth 
give rise to a ft force Tt which would then drive the computed orbit away from the true one. One can 
then use the variation of elements technique as given by Plummer to describe the effect of this 
force on the orbital elements and consequently the Cartesian coordinates.* 

A more direct and less complicated approach is the following. Denote the spherical coordi- 
nates of the satellite as r the radius, e the colatitude, and 4> the geocentric right ascension. Now 
the potential energy of the system is given by the standard solution of Laplace 1 s equation, or 
V = V(r,<9,0). The kinetic energy can be so written as 

T = +r 2 <9 2 +r 2 sin 2 6<p 2 ) . 

If the Lagrangian L is not an explicit function of time, then 


where the Hamiltonian 


• <9L 

H " / , q i dq k ~ L * 

i 

the sum taken over all canonical coordinates. Thus the Hamiltonian is a constant of the motion. 

In addition, if the forces present are derivable from a conservative potential, then H = E, the total 
energy. Assuming that the potential is a function of r and d, then, since the coordinate is cyclic, 
its canonically conjugate momenta = constant, or 

= mr 2 sin 2 Q 4> • 


*H. C. Plummer, "An Introductory Treatise on Dynamical Astronomy.’** New York: Dover Publications, 1960. 
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The Hamiltonian is 


H - E - ^mr 2 +"2 _ mr 2 ^ 2 +-2'mr 2 sin 2 ^(5{) 2 +V(r,^) 


- . ~2 mr 2 + ‘ V’ ( r, 0) , 


where ‘V’ (r ,9) is an effective potential. 

For stable motion along the radial direction, 

d ‘ V* ( r , 9) _ • (9V(r,9) V 

<? r " ' mrt + dr mr 3 sin 2 (9 


and 


a 2 ‘V ’ (r,t?) 
<9 2 r 


mti 2 + 


<3 2 V( r,t>) 

a / 2 


3Pj , 2 _ 

mr 4 sin 2 t. f 


> 0 , 


i.e. , this is a positive quantity, since V(r,^) goes as inverse powers of r, and therefore a stable 
point in the radius vector given by r 0 exists. Denoting 

d 2 ‘V’ (r,fc>) 


then, for a small perturbation, the radius vector will oscillate about the equilibrium position 
with a frequency given by 




hi | 1 0 2 V(r .g) 
, m dr 2 


op 2 


i 1/2 


i 2 r 4 sin 2 (9 , 


For the colatitude, let us assume a small change in angle ‘9 only so that 


2 ~ mr 2 9 2 + 1 V* ( r, 9) 


where 


‘ V’ ( r , (9) 


V(r,9) + 


V 

2mr 2 s in 2 5 
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and 


E - ‘ V * ( r , 0 ) 


If there were no displacement in 6, the relation between and the true colatitude, say 
9 - 6 q , would be given by 


d 1 V’ ( r , 9) 

d'9 


3V( r,6)\ *V cos 6 o 


d® Je Q mr 2 sin 3 6 Q 


Substituting P,, 


- 2 sin d 0 cos 0 0 \ d ~& 


dV( r,6)\ 


Then for u = 0, the expression for the total energy before displacement becomes 


sln "o /(9V(r, ) 

v < r ^o) 4 2 cosT7; l at— 


For an energy slightly larger than E 0 , and angular momentum given by P^, the angle ^ will 
oscillate about the value 0 O . Now setting 


S 2 ‘V’ (r,d) 
<9^ 2 


a 2 V(r.tf)\ (sin 2 e 0 * 3 cos 2 >\) / aV( r , u) 
d>< 2 ), + sin<5 0 cos«V \ d>' 


for small values of the displacement 0 - u 0 , ‘V’ (r,t>) can be expanded in a Taylor series such that 


‘ V ’ ( r , 0 ) = E 0 + -J k (d-^ 0 ) 2 , 


and the expression for the energy becomes 


E - E n 


= \ mr 2 0 2 + '2"k [d ~0 O ) 2 


This is in the form of the energy for a harmonic oscillator with coordinate ( 0 -& 0 ), mass mr 2 , 
and spring constant k . From this, the frequency of oscillation in 9 is given by 


JL Id 2 V (r,6>) 
ir 2 \ 96 2 


l_ ( sin2 ^0 + 3 co s 26 l o) / a v(r,g 
nr 2 sine o co S e o \ d& 


i 1 / 2 
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